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Abstract 

A general moment bound for sums of products of Gaussian vector's functions extending the moment 
bound in Taqqu (1977, Lemma 4.5) is established. A general central limit theorem for triangular arrays 
of nonlinear functionals of multidimensional non-stationary Gaussian sequences is proved. This theorem 
extends the previous results of Breuer and Major (1981), Arcones (1994) and others. A Berry-Esseen-type 
bound in the above-mentioned central limit theorem is derived following Nourdin, Peccati and Podolskij 
(2011). Two applications of the above results are discussed. The first one refers to the asymptotic 
behavior of a roughness statistic for continuous-time Gaussian processes and the second one is a central 
limit theorem satisfied by long memory locally stationary process. 

Keywords: Central limit theorem for triangular arrays; Moment bound for Gaussian vector's functions; 
Hermitian decomposition; Diagram formula; Berry-Esseen bounds; Long memory processes; Locally 
stationary process. 



1 Introduction 

This paper is devoted to the proof of two new results concerning functions of Gaussian vectors. The first 
one (Lemma [1] of Section [2]) is a moment bound for "off-diagonal" sums of products of functions of Gaussian 
vectors in a general frame. It is an extension of an important lemma by Taqqu (1977, Lemma 4.5). This 
result is useful for obtaining almost sure convergence and tightness of Gaussian subordinated functionals 
and statistics, see Remark [1] below. The proof of Lemma [1] uses the Hermitian decomposition of function 
and the diagram formula. A related but different moment bound is proved in Soulier (2001, Corollary 2.1). 



The second result is a central limit theorem (CLT) for arrays of random variables that are functions of 
Gaussian vectors, see Theorem [T] for a precise statement. Theorem [T] generalizes and extends earlier results 
due to Breuer and Major (1983), Giraitis and Surgailis (1985) and Arcones (1994, Theorem 2) to the case 
of non-stationary triangular arrays of Gaussian vectors. Extensions of the Breuer-Major theorem were also 
obtained by Chambers and Slud (1989), Sanchez de Naranjo (1993) and Nourdin et al. (2011). Most of 
the above cited papers treat the case of a single stationary Gaussian sequence and a function independent 

'Supported by a grant (No. MIP-11155) from the Research Council of Lithuania 



1 



of n. Generalization to stationary or non-stationary triangular arrays is motivated by numerous statisti- 
cal applications. Some examples of these applications, with a particular emphasis on strongly dependent 
Gaussian processes, are: statistics of time series (see for instance Bardet et ai, 2008, Roueff and von Sachs, 
2010), kernel-type estimation of regression function (Guo and Koul, 2008), nonparametric estimation of the 
local Hurst function of a continuous-time process from a discrete grid i/n,l < i < n (Guyon and Leon, 
1989, Bardet and Surgailis 2011, 2012). Two particular applications (limit theorems for the Increment Ratio 
statistic of a Gaussian process admitting a tangent process and a CLT for functions of locally stationary 
Gaussian process) are discussed in Section [SJ 

Starting with the famous Lindeberg Theorem for independent random variables, numerous studies devoted 
to CLT for triangular arrays under various dependence conditions had appeared. The case of martingale 
dependence was extensively studied in Jacod and Shiryaev (1987). Rio (1995) discussed the case of strongly 
mixing sequences. Some of more recent papers devoted to this question are Coulon-Prieur and Doukhan 
(2000) (with a new weak dependence condition) and Dedecker and Merlevede (2002) (with a necessary and 
sufficient condition for stable convergence of normalized partial sums). The CLT for linear triangular arrays 
was discussed in detail in Peligrad and Utev (1997) for several forms of dependence conditions. 

The case of Gaussian subordinated variables (functions of Gaussian vectors) is rather exceptional among 
other dependence structures since it allows for very sharp conditions for CLT in terms of the decay rate of 
the covariance of Gaussian process and the Hermite rank of non-linear function. These conditions are close 
to being necessary and result in CLTs "in the vicinity" of non-central limit theorems, see Breuer and Major 
(1981), Arcones (1994), Dobrushin and Major (1979), Taqqu (1979). The proofs of the above-mentioned re- 
sults rely on specific Gaussian techniques such as the Hermite expansion and the diagram formula; however, 
the recent paper Nourdin et al. (2011) uses a different approach based on Malliavin's calculus and Stein's 
method, yielding also convergence rates in the CLT. The main difference between our Theorem [1] and the 
corresponding results in Arcones (1994) and Nourdin et al. (2011) is that, contrary to these papers, we do not 
assume stationarity of the underlying Gaussian sequence (Y'„(fc)) and discuss the case of subordinated sums 
J2k=i fk.n{Yn(k)) whcrc /fc^n may depend on k and n. The last fact is important for statistical applications 
(see above). In the particular case when = / do not depend on k, n and (YVi(fc)) is a stationary process 
independent of n, Theorem [1] (iii) agrees with Arcones (1994) and Nourdin et al. (2011, Theorem 1.1). The 
proof of Theorem [T] uses the diagram method and cumulants as in Giraitis and Surgailis (1985). Section |4] 
obtains a Berry-Esseen bound in this CLT using the approach and results in Nourdin et al. (2011). Let us 
note that a CLT for Gaussian subordinated arrays is also proved in Soulier (2001, Theorem 3.1); however, 
it requires that Gaussian vectors are asymptotically independent and therefore his result is different from 
Theorems 

Notation. Everywhere below, X = {X^^\ . . . , X'^'^'>) designates a standardized Gaussian vector in M'', v > 1, 
with zero mean EX*^") — and covariances EX'^^^X'^"^ — Suv, u,v = I, . . . ,1^. Letter C stands for a constant 
whose precise value is unimportant and which may change from line to line. The weak convergence of 

T> 

distributions is denoted by — > . 

n— f 00 

2 A moment bound 

Let L^(X) denote the class of all measurable functions / = f{x),x = {x^^\ . . . , x'^'^^) € R"^ such that 
||/||2 := E/2(X) < 00. For any multiindex k = {k^^\ . . . ^k^"^) G Z"^ := {(j(^\ . . . , /'')) S Z'^,j(") > 
(1 < M < v)}, let Hf^{x) = (a;(^^) • • • iJfc(^) (a;^''^) be the (product) Hermite polynomial; Hk{x) := 
(— l)'^e^ /^(e~^ /2^(fc)^ /j — 0,1,... are standard Hermite polynomials (with (e^^ i]^q /j^j^ derivative 
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of the function x ^ e-^'/^). Write |fc| + . . . + fcM, fe! • • • fc^")!, k = (fc(i), . . . , fc''^') e Z^;. 

A function / e L-^(X) is said to have a Hermite rank to > if Jf{k) E/(X)iJ^(X) = for any 
fe e Z"^, |fe| < m, and Jf{k) 7^ for some fc, |fc| — m. It is weU-known that any / e L^(X) having a Hermite 
rank to > admits the Hermite expansion 

E ^^feW' (2-1) 



^ fc! 

|fc|>™. 



which converges in L^(X). 



Let (Xi,...,X„) be a cohection of standardized Gaussian vectors Xt = (Xj:^\ . . . , X^"^) £ W having 
a joint Gaussian distribution in M"^". Let e € [0, 1] be a fixed number. Following Taqqu (1977), we call 
(Xi, . . . , Xn) e— standard if \Exj:"^ xi^'' \ < e for any t ^ s,l < t, s < n and any 1 < u,v < v. 

As mentioned in the Introduction, Lemma [T] generalizes Taqqu (1977, Lemma 4.5) to the case of a vector- 
valued Gaussian family (Xi, . . . ,X„), taking values in R'^(;/ > 1). The lemma concerns the bound (|2.4p . 
below, where fi,t.m ■ ■ ■ , fp.t,n are square integrable functions among which the first < a < p functions 
fi,t,m ■ • ■ , fa.t,n for any 1 < t < n have a Hermite rank at least equal to to > 1 and where is the sum over 
all different indices 1 < ti < n [1 < i < p),ti ^ tj(i ^ j). In the case when fj^t.n — fj does not depend on 
t, n, the bound (|2.4p coincides with that of Taqqu (1977, Lemma 4.5) provided ma is even, but is worse than 
Taqqu's bound in the more delicate case when ma is odd. An advantage of our proof is its relative simplicity 
(we do not use the graph-theoretical argument as in Taqqu, 1977, but rather a simple Holder inequality). A 
different approach towards moment inequalities for functions in vector- valued Gaussian variables is discussed 
in Soulier (2001), leading to a different type of moment inequalities. 

Lemma 1 Let {Xi, . . . , X„) be a e— standard Gaussian vector, Xt = {X^^\ . . .,xi''^) e M'', > 1, and let 
.fj.t,n G L^(^), 1 < j < P > 2, 1 < t < n be some functions. For given integers to > 1,0 < a < p,n > 1, 
define 

g„ max ^ max lEX^^^^ X^-^r . (2.2) 

l<i<n ^ — ' l<u,v<v 
l<s<n,s^i 

Assume that the functions fi,t,m • ■ • , fa,t,n have a Hermite rank at least equal to m for any 7i > 1^1 < t < n, 
and that 

e<^—. (2.3) 

vp — 1 

Then 

'\E[fi.MAXt,)---fp,t,AXt,)]\ < Cie,p,m,a,i^)KnP-^Qt (2.4) 

where the constant C{e,p,m, a, v) depends on e^Pjmya,!^ only, and 

p 

l<t<n L .y, , 

Proof. Fix a collection (ti, . . . ,tp) of disjoint indices ti ^ tj{i ^ j), and write fj = fj^tj,m 1 < J < P for 
brevity. Let Jj{k) := Jfj{k) = E[/j(X)if^(X)] be the coefficients of the Hermite expansion of fj. Then, 

< ii/,iinEV2ij2^.,(x) 

i=l 

< ii/,iin(t'"!)"^=ii/,ii(*!)"^ 
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Following Taqqu (1977, p. 213, bottom, p. 214, top), we obtain 



|E/i(XtJ.../p(X,J| 



E E 



n 



< 



< 



1=0 \ki\+...+\kp\=2q U=i 

9=0 |fci| + ... + |fcp|=2g 



^^E E 

9=0 |fci|+... + |fep|=2g 



(fei!---fcp!)i/2 

,(ifc.i+...+ifc,i)/2En^^^^^ni<,-<,vw 

(fci!---fcp!)i/2 



< i^i^ ^ (e(i.p~l))(l'^il+-+l^-l)/2 <oo, 

9=0 |fci|+... + |fep|=2g 

where X - 7V(0, 1) and 

i^i ||/l,ti,n||---||/p,tp,«|| < i^, 

where X is defined in (j2.5l) and X is independent of ti, . . . , tp, and where we used the assumption 
get the convergence of the last series. Therefore, 

oo 



to 



E E' 



9 = Ifcil + . . . + Ifcpl = 2q 
fell > Ifcal > m 



|Egfc^(X,J...gfc^(X,J| 
(fei!---fcp!)i/2 



Now, the following bound remains to be proved: for any integers m > 1,0 < a < p,n > I and any 
multiindices fci, . . .kp G Z^j^ satisfying |fci| + . . . + \kp\ = 2q, |A;i| > to, . . . , \ka\ > to, 

J2'm,^S^,,)---H,^^{X,J < Ci(£(i^p-l))(l'^^l+-+l^-l)/2(fci!...fc,!)i/2„P-fQ|, (2.6) 

where Ci is some constant depending only on p, v, a, e, and independent of fci, . . . , kp, n. 

First, we write the expectation on the left hand side of (12.61) as a sum of contributions of diagrams. Let 



T 



^(1,1) (1,2) ... 
(2,1) (2,2) ... (I,fc2) 



(2.7) 



\{p,l) (p,2) ... {p,kp)J 

be a table having p rows ti, . . . ,Tp of respective lengths |t„| = fc„ = |A;„| = fci^-* + . . . + fcl'^^ (we write 
T — U^=i '^«)- A sub-table of T is a table T' — [J^f^u '^u, U C {1, . . . ,p} consisting of some rows of T written 
from top to bottom in the same order as rows in T; clearly any sub-table T' of T can be identified with a 
(nonempty) subset U d {1, . . . ,p}. A diagram is a partition 7 of the table T by pairs (called edges of the 
diagram) such that no pair belongs to the same row. A diagram 7 is called connected if the table T cannot 
be written as a union T ~ T' U T" of two disjoint sub-tables T', T" so that T' and T" are partitioned by 7 
separately. Write r(r),rc(T) for the class of all diagrams and the class of all connected diagrams over the 
table T, respectively. Let 

p(t,s):= max \EX[^'^ X^^\ (t^s). 

1<U.V<U 

Note < p{t,s) < e and Q„ = ''^^'^i<t<nJ2i<s<n s=^t P"^(^^^)- ^he diagram formula for moments of 
Hermite (Wick) polynomials (see e.g. Surgailis, 2000), 



\EH,^^iXt,)---Hf^^{XtJ < E n iPitu,t.)Y- 

7er(T) i<u<v<p 

= E n E n (pitMY 

(Ui,...,Uh) r=l leFaiUr) u,veUr,u<v 



(2.8) 
(2.9) 
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where is the number of edges between rows r„ and r„ in the diagram 7 over table T, and the sum 
'^{Ui Uh) taken over ah partitions {Ui, . . . , Uh), h — 1,2, ... , [p/2] of {1, . . . ,p} by nonempty subsets 11^ 
of cardinality \Ur\ > 2. (Thus, (|2.9p follows from (I2.8P by decomposing 7 e r(T) into connected components 
7r € Tc{Ur),r = 1, . . . ,h;h = 1, . . . , [p/2]; the restriction \Ur\ > 2 stems from the fact that any edge must 
necessarily connect different rows.) From (|2.9|) we obtain 

h 

5]'|E[/i(X,J.../,(X,J]| < E n E ^".^^^(t)' (2.10) 

(c/i,...,c/fe)r=i^er,(c/,) 

where, for any sub-table U C T having at least two rows and for any connected diagram 7 G Tc{U), the 
quantity In,uil) is defined by 

ii,vtE U,u<v 

where (recall) the product is taken over all ordered pairs of rows {tu,t^),u < v oi the table U, and luv is 
the number of edges in 7 between the wth and the wth rows. Below we prove the bound 

InAl) < if3^l'^^'l/^r^l^l-^(nQ„)^, (2.11) 

where \kij\ := J2u<£U ^he number of points of table U and a{U) :— |{1, . . . , a} P| [/| = #{u G U : |fc„| > 
to} is the number of rows in U having at least m points. Clearly, it suffices to show (j2.11l) for U ^ T. 
Next, let for 1 < u,v < p,u ^ v, denote 

^ X fc„//c„\^-/'=" 

El<,< Ei<< ^y-i^^t)) ) ■ (2-12) 

Let A :— {1, . . . , a}. A' :— \\, . . . = {a + 1, . . . It follows immediately from the definition of R^v 

and pis, t) that 



R11 



By the Holder inequality (see Giraitis and Surgailis, 1985, p. 202, for details), 



(2.13) 



In,T{l) < min 

n n • (2.14) 

For any subset [/ C { 1 , . . . , , let 

LiU):=Y. E If' ^*(t^)-E E If' (2-15) 

u^U u<v<.p u£U l<v<u 

L := L(T),L* := L*(r). Clearly, 

L([/) + L*([/) = E^ E ^™ = 1^1 (2.16) 

u£U v—l,...,p,v^u 

is the number of points in U. From (I2.13P - (|2.14l) . 

/n,T(7) < min (^i*+i(^^)Qi(^)£|T|/2-n^L(A)^ ^L+i*(A^)gL*(A)^|T|/2-™L*(A)^ 

where |r| = ^^^-^ fc,. As < L(A),L*(A) < p, see (EH]), we obtain 

In,Th) < el^l/'-™^min(n^'(^)+^*(^^)+^(^^)g^(^),n^(^)+^(^^)+^*(^^)Qf (-4)) 

= el^l/2-™PnP-"niin(n^*(^)g^(^),n^(^)gf f-^)) 

= £l^l/2-™Pn''-"(ng„)^ min ((n/g„)^"^^^^ , (n/Qnf^^^^' 
< el'^l/2-"PnP-"(ng„)t, 
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proving (|2.1ip . 

With (P1T|) - (PAU)) in mind, 

5:'|Ei7,^(X.J...Ffe^(XJ| < C3sl-I/^ ERE m^^^'^Q^ 

{Ui,....,Uh)r=i',er^(Ur-} 

= Ce^^y^n^-^Ql E n E 1 
{Uu---:Uh)r=i-fer,{Ur-) 

= Cael^l/^n^'-tQl ^ 1, 
7er(T) 

where the last sum {— the number of aU diagrams over the table T) does not exceed 

lEiJ^a) (X)--- (X) . . . (^) • • • {X)\ < {pv - l)i\k.\+■■■+\K\)n^k^\ ■ ■ ■ k^)^/\ 
see Taqqu (1977, Lemma 3.1). This proves the bound ()2.6p and the lemma, too. □ 

Lemma [1] can be extended to non-standardized Gaussians as follows. To this end, we introduce some 
definitions. Let Y = {Y^^\ . . . ,Y^'^'>) S be a Gaussian vector with zero mean and non-degenerate 
covariance matrix S = (EY'-"^Y'^^^)i<u.v<u- Let L^(l^) denote the class of all measurable functions / : 
R'' — > M with Ef'^{Y) < oo. Let m > be an integer. We say that / e L,'^{Y) has a generalized Hermite 
rank not less than m if either m = 0, or to > 1 and 

E[P{Y)f{Y)] = for aU P e V,n-i{^l (2-17) 
hold, where Vm{^'^) stands for the class of all polynomials P in variables y'-'^', . . . , j/'"^^ of degree m, that is, 

Piy) = Eo<iji<™ <j)y^ = E,<i)>o....,,(^)>o:,(i)+...+,<.^)<™ • • • - j^'^') (y^'^P''' ■ • • (y^^^V'"" • 

Let X := E-i/2y, f{x) := f{Y^I'^x). Then X has a standard Gaussian distribution in W and / e ^■{X) 
with 

||/||2 = E|/(X)|2=E|/(y)p. (2.18) 

The following proposition is known, see Nourdin et al. (2011, Proposition 2.1), Soulier (2001, p. 195), but 
we include a proof of it for completeness. 

Proposition 1 Let Y^X, f G L^(l^), / G L^(X) be defined as above and m > he a given integer, f has 
a generalized Hermite rank not less than to if and only if f has a Hermite rank not less than m. 

Proof. The above proposition is true ii Y — X has a standard Gaussian distribution; see Soulier (2001, 
p. 194). By definition 

E[PiY)f{Y)] = E[P(X)/(X)], (2.19) 

where P{x) P{j:'^/'^x). Clearly, P e Vm-ii^") implies that P e 7',„-i(M'") is a polynomial of degree 
TO — 1. Therefore / having a Hermite rank not less than to implies by p.l9p that / has a generalized Hermite 
rank not less than m. The converse statement again follows from ()2.19p . by taking P{y) = P(E~^/^y), where 
P G Vm-ii^") is an arbitrary polynomial of degree to — 1. □ 

Let (1^1, . . . , Yn) be a collection of Gaussian vectors Yt — {y}^\ . . . , !"/"•') e with zero mean EYt — 
and non-degenerated covariance matrices St = (Cov(y/"\ ^/"■'))]^<„ „<j,j having a joint Gaussian distribu- 
tion in W^. Let e e [0, 1] be a fixed number. Gall (Yi, . . . , y„) e-correlated if |Cor(r/"\ ri"^) | < e for any 
t s,l < t, s < n and any 1 < u,v < ly. Clearly, if the Yt's are standard, this is equivalent to (l^i, . . . , Yn) 
being e— standard. 
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We also use some elementary facts about matrix norms. Let \x\ = {J2i=ii^''^^)^Y^^ denote the Euclidean 
norm in M"^, A = (aij) a real x i/— matrix, the transposed matrix, / the unit matrix, and ||A|| := 
^uP|cc|=i l^^l the matrix spectral norm, respectively. Then ||A||oo '■— T^sxi<i,j<u \o-ij\ < ||^|| < ^^H^Hoo and 
< ||A|j ||i?|| for any such matrices A,B. An orthogonal matrix O = (oij) satisfies OC' = O'^O = I 
and ||0|| = = 1. Any symmetric matrix A can be written as A — O'^AO, where O is an orthogonal 

matrix and A is a diagonal matrix. In addition, if A is positive definite, then = j|Aj| = Amax, ll^^^ll = 
||A^^|| = A~[jj, where A,nax > ^mm > are the largest and the smallest eigenvalues of A. We shall also use 
the facts that for any symmetric positive definite matrix A, 

||Al/2|| = j|A||V2^ P-1/2|| ^ p-l||l/2^ (2.20) 



since 



|Ai/2|| ^ II0TA1/20II = ||Ai/2|| = xUL = IIOTA-V20II = ||A-i/2|| = x-^/\ 



Corollary 1 Let {Yi, . . . ,Yn) be an e— correlated Gaussian vector, Yt ~ {Y^^^\ . . . ^Y^'^^) £ {i/ > 1), 
with zero mean FiYt = and non- degenerated covariance matrices St satisfying 

max < Cniax (2.21) 

l<t<n 

for some constant Cmax > 0. Let fj,t,n € ^^(l^j), 1 < j < p {p > 2),1 < t < n be some functions. For 
given integers m > 1, < a < p, n > 1, let Qn denote the sum in (dl) where X^''\ xi"^ are replaced by 
y}'^\Ys^\ respectively. Assume that the fund 
equal to m for any n > 1,1 < t < n, and that 



y}'^\Ys^\ respectively. Assume that the functions /i,t,n, • • ■ ,fa,t,n have a generalized Hermite rank at least 



e < . (2.22) 

(vp- l)z^^Cmax 

Then 

Y!mfi,tun{yt,)---fp.t,.n{'^t,)]\ < CKnP-^Ql, 

where K := 11^=1 niaxi<t<„ E^/^ |^j2^ ^^CKj)] and the constant C = C{e,p, m, a, v, Cmax) depends on e,p, m, a, 
i^, Cmax only. 

Proof. We will reduce the above inequality to that of Lemma[Tl as follows. Let Xt := ' Yt, fj.t,n{x) := 
fj^t.ni^y^^)- The Xt's have a standard Gaussian distribution in R'' and the fj^t.nS, satisfy = 
E[/2j ^(l^t)] , see (|2.18l) . By Proposition [1] fj^t.mj = 1, ■ . ■ ,a have a Hermite rank not less than m. Next, 
using (|2.20p . (|2.2ip and the fact that the Yt's are e— correlated, for any t ^ s, 1 <t, s < n,l < u,v < v 

|EA(")a(-)| < i.2||s-i/2|,^|,5.-i/2|,^ |Ey/"VW| < eA\^-;"\\^-^l^\\ < e^^,,^ (2.23) 

l<u,v<iy 

This implies that the Gaussian vector {Xi, . . . ,X„) e M"" is e— standard, where e :— ei^^Cmax- Then, in 
view of (|2.22[) . (|2.4p of Lemma [T] applies, according to which 

|E[/l,ti,„(FtJ • • • /j,,tj,,„(Ytj,)] I = ^ \E[fi^ti.niXti) ■ ■ ■ fp,tp,niXtp)]\ 

< C{e,p,m,a,v)KnP~'^Qi < C{£,p,m,a,v,Ci^s.^)Kn^~^ Qj , 

where K,Qn are the corresponding quantities in Lemma [T] (|2.4p satisfying K = K,Qn < Cmax)"^ Q n by 
(PrSj) . (11231), respectively. □ 

We remark that condition (|2.22p is not optimal since it does not reduce to (|2.3p in the e— standard case. 
This loss of optimality is due to the use of robust inequalities for matrix norms in (|2.23p . 
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Remark 1 As mentioned in the Introduction, Lemnia[T]and Corollary[T]can be used for proving the tightness 
and the strong law of large numbers of various non-linear statistics from Gaussian observations. See Bardet 
and Surgailis (2011, 2012) on application for roughness estimation and Csorgo and Mielnichuk (1996), Koul 
and Surgailis (2002) for empirical process. The above-mentioned applications concern the 4th moment bound 
Mn := E(^"^j ft,n{Ynit)))^ = 0{n~'^) for a suitable k > 0, where (l^„(t)), (/t,„) satisfy similar conditions 
as in CoroUary [TJ Clearly, M„ = 44=1 E[ 11^=1 /ti,"(^n (*«))] can be decomposed into four terms 

according to the number of coinciding "diagonals" ti — tj in the last sum, where each term can be estimated 
with the help of Corollary [T] Let us note that condition (I2.22p in the above applications is guaranteed by a 
preliminary "decimation" of the sum /t.n(^n(i), see (Csorgo and Mielnichuk, 1996) and (Bardet and 

Surgailis, 2012) for details. 



3 A CLT for triangular array of functions of Gaussian vectors 

Let (X„(fc)) 

i<fe<n.neN ^ triangular array of standardized Gaussian vectors with values in M'', X„(fc) — 
ixi^\k), . . . ,xi"^(fc)), ExiP\k) = 0, ^X^,r\k)xi'\k) = 5„. Now define, 

r(f''^)(j, fc) EX(P)(j)X(«)(fc) (1 < J, fc < n). 

For a given integer m > 1, introduce the following assumptions: for any I < p,q < i'-, 

sup max ^ |r^«)(j,fc)r < 00, (3.1) 

1 I m 

sup-^ riP-«)(j,fc) ^0, (3.2) 



n>i ri 

— I < j, k < n 

\j - k\ > K 

V(j, fc)e{l,...,n}2, |r(f'')(j,fc)| < \p{j-k)\ with Y.\p{j)r<^. (3.3) 

jei. 

Note dSSl) (EH) and 1^ =^ Let L2(X) := {/ e lJ^{X) : Ef{X) = 0}, where X denotes a 

standard Gaussian vector as above. 

Theorem 1 Let (-X^n(fc))i<fc<„ be a triangular array of standardized Gaussian vectors. 

(i) Assume h3.1\) . Let fk € L§(JC) (1 < fc < n) have a Hermite rank at least m £ N* . Then there exists a 
constant C independent of n and fkj^ l£ k < n such that 

e(""'/'E^-(^"W)) < max (3.4) 

\ ^ — ' / l<k<n 

fe=l 

(ii) Assume i3. 1]) and \3.2]) . Let fk^n G 1Lq(X) (n > 1, 1 < fc < n) he a triangular array of functions all 
having Hermite rank at least to € N*. Assume that there exists a \J^{X)^valued continuous function 
(/)t,t € [0, 1], such that 

sup ll/[r«],« - = sup E(/[^„]^„(X) - 0^(X))^ — > 0. (3.5) 

Te(04] Tg(0,l] n-foo 

Moreover, let 

al E(n-i/2 ^ /, „ (x„(fc)) ) a\ (3.6) 



n^oo 



where a > 0. Then 

n 

n-^'^y^h,n{Xn{k)) ^ AA(0,a2). (3.7) 
fe=i 
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(iii) Assume \3. Moreover, assume that for any t e [0, 1] and any J €z N* , 

{Xni[nr]+j))^ A {WAj))_j^^^j, (3.8) 

— '— n— !-oo — ' — 

where iWr{j))j^x stationary Gaussian process taking values in and depending on parameter 
T G (0,1). Let fk.n € Lq(X) {n > 1,1 < k < n) satisfy the same conditions as in part (ii), with 
exception of fi3.6\) . Then LS. 6]) and Jg. 7p hold, with 



(3.9) 



We remark that parts (i) and (ii) of Theorem [T] are natural extensions of Theorem 2 of Arcones (1994) 
(for instance, condition (13. ip is the same as condition (2.40) of Arcones (1994) in the case of stationary 
sequences). We expect that parts (i) and (ii) can be also obtained following the method in Nourdin et al. 
(2010). Part (iii) seems more interesting. In Bardet and Surgailis (2011), (iii) is applied when Xn{i) ~ ^jjn 
and ^Zt)t is a vector valued continuous time process. 

Similarly to Lemma [H Theorem [T] can be extended to nonstandardized Gaussian vectors. Corollary [2] 
below refers to the most interesting part (iii) of Theorem [1] 

Corollary 2 Let Yn{k) = (vj-^^k), . . . ,Yj:''\k)) e W , 1 < k < n, n e N be a triangular array of 
jointly Gaussian vectors, with zero mean EY'„(fc) = and non-degenerate covariance matrices Sfc.„ = 
YiYn{k)Yn{ky . Assume that covariances rn'''\j, k) := Cov(yi^^(j), ri«^(fc)) satisfy lOI) . for some m>l. 
Moreover, assume that LS. 8\) holds with Xn{-) replaced by Y„{-), ivhere (WV(j))j-g2 ^ stationary Gaussian 
-valued process with non- degenerate covariance matrix Et- :— EWV(0)WV(0)''' such that 

sup ||E;1|| < oo (3.10) 
Te(o,i] 



and 



sup I IE 1 1 — > 0. (3.11) 

Te(0,l] n-J-oo 



Let fk.n € lLQ(Y'„(fc)), l<k<n, nENhea triangular array of functions all having a generalized Hermite 
rank not less than m and such that 

sup e(/V„],„(X)-0,(X))' 0, (3.12) 

Te(0,l] ^ ' n-J-oo 

where fk.n{x) :~ fk.n(^]/nX) and where 't>T,T G [0,1] is a lJ^{X) — valued continuous function, with X a 
standard Gaussian vector in M.^ as usual. Then 

n 

n-'/'y2fk,n{Yn{k)) A AA(0,a2). (3.13) 



fe=i 



where is defined in i3.9\) . with (I)t{x) :— 4>t{^t^^'^x). 

1/2 

Proof of Corollary\^ Similarly as in the proof of CoroUaryll] let X„(fc) :— E^, J y„(fc). The JC„(fc)'s are 
standardized Gaussian vectors in W and the fk,nS belong to 'L?{X) and have a Hermite rank not less than 
m. Assumptions p.8p and p. lip entail for any r e (0, 1), (X„(j + [nT])))_ {W^[j)^ _ 

— —1/2 

where WV(j) '■— Er W,-(j), j G Z is a stationary Gaussian process having a unit covariance matrix 
¥Wr{Q)Wr{Qy = I. Conditions (jXTUI) and (jXTTI) imply that maxi<fe<„ ||E^y^|| < C. The last fact to- 
gether with condition p.3p for covariances r^''^\j, k) := Go\{Yn^\j), Yn''\k)) imply a similar condition for 
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Gov{X^''\j),X^''\k)): for all {j, k) e {1, • • • , n}^ we have that maxi<„,„<, \&xt\i)x},^\k)\ < C\p{j~k)\; 
see (|2.23l) . This way we see that the conditions of Theorem [IJiii) including l|3.5p are satisfied and can be 
applied to the families of Gaussian vectors (X„(fc)) and functions [fk,n), yielding p.l3p . □ 

Proof of Theorem[l[ (i) Using Arcones' inequality (see Arcones, 1994, (2.44) or Soulier, 2001, (2.4)), one 
obtains 



n 2 1 " 1 / 

e(""'/' E (^'-(k)) ) = - E ii/fcii' + - E EA(x„(fc))/,(x„(£)) 



n ^ — ' n 

fe=l k=l 



< max ||/fe||2 + C( max ||/fc||) max V max \AP-'^^k,e)\ 

- l<k<n"-' " M<fe<n l<k<n ^ \<v,q<v^ " ^ ' 

^ ^ ^ ^ ^ ^ \<.i<n,l=^k 



where C is a positive real number not depending on n or fk- Now, using assumption (jSHJ, (i) is proved. 

(ii) We use the following well-known fact. Let (Z„)„>i be a sequence of r.v.'s with zero mean and finite 
variance. Then — >■ A/'(0,cr^) if and only if for any e > one can find an integer no(e) > 1 and a 



sequence {Zn^t)n>\ satisfying Zn,t — > A/'(0, Cg) and Vn > no(e), E(Z„ — ^n,c)^ < £• 

n— >C30 

Let Zn :— n-^^^^ '^k=i ff^,^ i^nik)). We shall construct an approximating sequence Zn,e with the above 
properties in two steps. 

Firstly, by condition p.Sp and continuity of (pr, for a given e > one can find integers M,no{e) and a 
partition —: tq < ti < . . . < tm < tm+i '■= 1 such that V n > no(e), 

1/2 

max max ||/fc,n - 0tJ| = max max (E(/fc „(X) - 0.^, (-^))^) < e- (3-14) 

Put 

M 

Z„,, :=n-i/2^ E (X„(fc)). 

i=0 fc/Tie(ri,ri+i] 

Note for any r £ (0, 1], the function V'r has Hermite rank not less than m, being the limit of a sequence of 
Lq(X)— valued functions of Hermite rank > m. Therefore for the difference Z„ — Zn.e the inequality p.4p 
applies, yielding M n > no{e) 

E(Z„-Z„.,)2 < C max max ||/fc,n - 0rj|' < Ce' (3.15) 

0<i<M k/nG{Ti,Ti^i] 

in view of p.l4p . with a constant C independent of n, e. 
Secondly, we expand each 0t-. in Hermite polynomials: 

^n{x) = E '^^ki^^' (z = 0,l,...,M) (3.16) 

m<|fc| 

where 

J,(fc) J^^^(fc) =E0,.(X)ff;^(X), |J,(fc)| < 110,, II (fc!) 1/2. 
We can choose t(e) G N large enough so that 

||0r. -0n,e|| <e, (i = 0,1,...,M), (3.17) 

where 4>Ti,e is a finite sum of Hermite polynomials: 

cPr^x) := E (* = 0,1,...,M). (3.18) 

Tri<|fc|<t(e) 
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Note i(e) does not depend on i = 0, 1, . . . , M, and e > is the same as in p.l4p . Put 

M 

Z„,e:=n-i/2^ 0n,e(X„(fc)). (3.19) 

i=0 fc/ne(Ti,Ti+i] 

Applying p.4p to the difference Zn,e — ^n,e and using p. 171) and (I3.15p . we obtain V n > no(e), 

E(Z„-Z„,,)2 < Ce" (3.20) 

where the constant C is independent of n, e. Let ^ := EZ^ ^. From p.20p and condition p.6p it follows 
that Vn > nQ{e), 

a'^-Ce < cr^ , < cr^ + Ce, (3.21) 

with some C independent of n, e. In particular, by choosing e > small enough, it follows that lim inf„_>.oo ^ > 
0. We shall prove below that for any fixed e > 0, 

7 1 

t/„ — = E 'l^r.AXnik)) A AA(0,1). (3.22) 

As noted in the beginning of the proof of the theorem, the CLT in ((X7|) follows from (15^ . (jX^ . ((X^ . 
Indeed, write 



Ee''^^" - e-''''"'/^ = (Ee'"^" - Ee'''^"-^) + (Ee'""^-''^" - fT'^^^J'^ 

i=l 

Here, for some constant C independent of n, a, e, 

\eiin)\ < Ei/2|c-(^"-^-=)-lf < |a|Ei/2|z„-Z„.,|2 < qa|e, 

and therefore £i{n),i — 1,3 can be made arbitrarily small by choosing e > small enough; see (I3.20p . p.2ip . 
On the other hand, the convergence in p.22p implies uniform convergence of characteristic functions on 



for any A > 0. This 



compact intervals and therefore sup|^|<^ K2("-)| < sup|^|<2yi 
proves p.7p . 

It remains to prove (j3.22p . The proof of the corresponding CLTs for sums of Hcrmite polynomials in Ar- 
cones (1994) and Breuer and Major (1983) refer to stationary processes and use Fourier methods. Therefore 
we present an independent proof of (|3.22l) based on cumulants and the Holder inequality in (I2.14p . Again, 
our proof appears to be much simpler than computations in the above mentioned papers. 

Accordingly, it suffices to show that cumulants of order p > 3 of f/„ asymptotically vanish. In view 
of (|3.2ip and linearity of cumulants, this follows from the fact that for any p > 3 and any multiindices 
fe„ = {ki^\ . . . , ki"^) e Z";, u=l,...,p with fc„ = |fc„| = /ci^' + . . . + ki;"^ >m {l<u<p), 

n 

S„ := Yl |cum(ii,...,tp)| = o(nP/2), (3.23) 
ti,...,tp=i 

where cum(ti, . . . ,tp) stands for joint cumulant: 

i(ii,...,ip) cum(iJfe^(X„(ii)),...,fffc^(X„(tp))) . (3.24) 
Split S„ = S'„(if) + S;^(i^), where 

n 

ti,...^tp—i 



"PI 

cum I 
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and where K will be chosen large enough. Then for any fixed K , we have T.'^{K) = 0{n) = o{vF/^) as p > 3. 
The remaining sum does not exceed '}2i<i=^j<p^n,i,j{^)^ where 

n 

K,^AK) ■■= E \cnm{h,...,tp)\l{\U-t,\> K). 
ti ,...,tp— 1 

Therefore, relation p.23p follows if we show that there exist S{K) — > and no such that for any 1 < i 

K—>oc 

j < p and any n > fiQ 

limsupl];; ,j(i^) < S{K)nP/^. (3.25) 

n— fcjo 

The proof below is limited to = (1, 2) as the general case is analogous. It is well-known that the joint 
cumulant in (I3.24[) . similarly to the joint moment in (12. 6p . can be expressed as a sum over all connected 
diagrams 7 G Tc{T) over the table T in (|2.7p . By introducing p{s,t) := maxi<p^g<i, rn''^\s,t) , we obtain 



|cum(ti,...,tp)| < 51 n {p{iu,t^)Y-, (3.26) 
7erc(T) i<u<v<p 



where we use the notation in (j2.6p . Therefore, 

n 

KiAK) < E E n {p{tu,t,)Y-H\ti-t2\>K) E ^".^(7)' 

7er<,(T) ti,...,tp=i i<u<v<p 7erc(T) 
Next, by applying the Holder inequality as in ()2.14p . 



In,T{l) < min 

.l<u<v<p l<u<v<p 



where (cf. ^J2\i ) 



{s,t)f^"^^ (^,z;)^ (1,2), (2,1), 

Ei<*<„(Ei<.<„p'^(s,i)l(|i-s|>^)) j , {u,v)^{ia), 
Ei<t<„(Ei<.<„p'^(i,s)l(|i-s| >^)) j , (w,z;) = (2,l). 

From assumptions (13.11) . (13. 2p , there exists a constant C and — > independent of n such that for 
any k > m and any n > 1 

n 

sup y^p^{s,t) < Cn, 

l<t<n . 
— — s—\ 

n 

sn^ ^p^{s,t)l{\t- s\> K) < 5{K)n. 
1=1 

(Cn'.Jk.^ (u,«)^(l,2),(2,l), 
Ruv < < J(if)n^"/'=^ {u,v) = (1,2), 
[ (5(ii:)n^i2/'''i, (u,w) = (2,l), 

with some 5{K) — > independent of n. Consequently, the minimum on the right-hand side of (13.27^ 
does not exceed 

C~5{K) min (n^^<^<^<. '-^'^^ , n^^<^<^<^ = C5(K)n™"(^(^)'^*(^» 



Therefore 
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where the quantities L{T),L*{T) introduced in ^^^\^ satisfy L{T) + L*{T) = p, sec (^l^ . and therefore 
niin(L(r),i*(r)) < p/2. This proves ((X^ and the CLT in (jX^ . thereby completing the proof of part 

(ii) . 

(iii) Let us first prove p.6p with cr^ given in p.9p in the case when fk^n = / do not depend on fc, n (in such 
case, one has 1/)^ = /, too). We have 

al = n-^ y- El/(X„(fc))/(X„(fe'))| = / ^^n^dr, (3.28) 



n „i 

-1 ^ E[/(X„(fc))/(X„(fe'))] = / 



where 

n— [nr] 

^»(t) ^ E[/(X„([nr]))/(X„(M+j))]. (3.29) 

j = l-["T] 

Condition (13. 8p imphes that 

E[/ (X„([nr])) / (X„(M + j)) ] ^ E[/ (W.(0)) / (Ty.(j-)) ] 
for each j G Z as n — )■ 00. From (13. 3p and with the inequahty of previous part (i), it exists C > such that 

E[/(X„([nr]))/(X„(M+j))]| < C\p{j)r , (3.30) 
and X^jgz |p(j)r" < Hence, from Lebesgue Theorem, 

K(r) =^l,e{i-[„.],..,„-[„.]}E[/(X„([nr]))/(X„([nr]+j))] ^ E[/ (T^.(O)) / (W,(j-)) ] ■ 

The dominated convergence theorem ahows one to pass to the hmit under the integral, thereby proving p.6p 
with (T^ given in p. 91) in the case /fc^„ = /. 

To end the proof, consider the general case of fk,n as in (iii). Let Zn.t be defined as in p.l9p . Note relation 
p.20|) holds as its proof does not use (|3.6p . In part (ii), we used p.6p to prove (I3.2ip . Now we want to prove 
p.2ip using p.Sp instead of p.6p . This will suffice for the proof of (iii), as the remaining argument is the 
same as in part (ii). 

Consider the variance cr^^j = EZ^ ^ of Z^^t defined in p.l9p : 



-1 I ^ EA' + 2 EAi^, 

iO<i<M 0<i<j<M 



where 

k/ne[Ti,Ti + i) 

Let us show that for e, Af fixed, and as n — > 00, 



EAA = o{n) (i^j), (3.31) 

n-iEA' ^ / VE[0,,.,(W"r(O)).^.,.e(Wr(j))]dT. (3.32) 

Here, p.32p follows from the argument in the beginning of the proof of (iii), as 4'Ti,e does not depend on 
k, n. Relation p.3ip is implied by the following computations. Using the Hermitian rank of functions (j)Ti,e, 
for i < j one obtains 

|E0,^,,(X„([nT,]+fc))0,^..,(X„([nr,]+^))| < C^r^A ■ Ur,.c\\ max r^f ([nr.] + fc, [nr,] + £) 

< C\\4)r,J ■ Ur,d\p{[nTj] ~ [nr,] + ^-fc)|'". 
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Therefore, for i < j, and e small enough, 



[Ti+in]-[Tin] [Tj + in]~-[Tjn] 



\ED,Dj\ < C max \\(j)r\\^ ^ ^ \p{[nTj] - [nn] + 1 - k)]"" 

n 

< C max lU^lP Vfc|p(fc)r o(n) 

re[0.1] II ^ 'I 
/c— 1 

since ELi^|p(^)|" ^ \/^Ei<fc<v7r + "Efc>vH |p(^)r = ''H- Thus, p.3ip is proved. From 

([33T|1 . (|332l) it follows that for any e > 

lim = := ^ / ^ E [0..,, (VF,(0)) 0..., (PF.(j)) ] dr. 



Consider the difference a'^ — a"^ — X]i=o /J^'^^ 5I)jez QM,e(''', j)dT, where 

< + (3.33) 

Using uniform continuity of 0r, t G [0, 1] (in the sense of L^— norm continuity), we obtain that the right-hand 
side of p.33p can be made arbitrarily small by choosing M {— the number of partition intervals of [0,1]) 
and t{e) {— the truncation level of Hermite expansion) sufficiently large, uniformly in r G [0, 1] and j G Z. 
On the other hand, \&M,€{T,j)\ < C sup^f^^Q ^ ll0T|P|p(j)|™ by Arcones' inequality, c.f. p.30p . Therefore 
\^M,eiT, j)\ is dominated by a summable function uniformly in M, e. Now, p.21[) follows by an application 
of Lebesgue theorem. This proves part (iii) and Theorem [T] too. □ 

4 A Berry-Esseen-type bound for nonstationary Gaussian subor- 
dinated triangular arrays 

This section obtains a Berry-Esseen-type upper bound in the CLT (13. 7p for non-stationary Gaussian subor- 
dinated triangular arrays following the method and results presented in Nourdin et al. (2011). We will refer 
NPP to the last paper in the rest of this section. To simplify the discussion, we restrict our task to the case 
when the functions fk,7i = / in Theorem [T] (iii) do not depend on k,n. As in NPP, our starting point is the 
Hermite expansion (|2.ip written as 

oo 

/ = E /W' E Jf(^)Hklk'- (4-1) 

'^=™ \k\=i 

Following NPP and using the Hermite expansion in (14. ip . we first define the following quantities: for j £ Z, 
l>m,N>m,n£W and J G {1, . . . , n}: 

0{3) ■■= \pU)\, K:=M{keN:dU)<-,m>k}, 0:=V0(j)", (4.2) 

n 

^In ■■= ^i"' E Cov(/(,)(X„(t)),/(,)(X„(t'))), (4.3) 

t.t' — — 71 

- (2^ E ^(^T E (^(^y^") ^'^ (for 1 < e < ^ - 1), (4.4) 



b1<" li'l<" 



1/2 

^2,Ar 2(2K^v^-B){¥.\f{X)\ E E[/2)(X)]) , (4.5) 



14 



N , _ p e-1 2 



^3,„,^ := Infix)] ^ E J'J- ( ■ ) v/(2^ - 2j)!7«,o) , (4-6) 



^6.„,J := Jia/IL sup ^ ||E[X„(M)^,I(M+j)]-E[Ty.(0)W-^(j)] , (4.9) 

^ 0<T<1. Jl 



|J|<./ 



A 



7, J 



iE[/2(X)]i.™ ^ 0™(fc). (4.10) 



|fe|>J 



Note that terms v42,n, ^3.ri,Af and ^5,„.Ar are the same as m NPP, A4,n,N is a minor improvement of the 
corresponding term in NPP, and Ag.n^j reflects the "convergence rate" in p.Sp . Term yli.„ of NPP (which 
does not appear in our bounds) is "absorbed" in the term inf i<j<„ ^7^,/ in the bounds (i)-(iii), below, due 
to a somewhat a different approximation (see (I4.15p ). 

Proposition 2 Let the assumptions of Theorem{I\ (in) prevail, with fk^n = f for all 1 < k < n, n (z N* 
where f : R is a Lipschitz function with \f{x) — f{y)\ < |i9/|oo \x — y\ for all x ^ y (z . Define 

S„ n-i/2^"^i/(X„(t)) and let S be a zero-mean Gaussian random variable with a variance ai := 
/o CoYifiWrmJiWrU))) dr < c^. Then: 

(i) For any function h twice continuously differentiable with bounded second derivative and for every 
n> K, 

E[h{Sn)] -E[h{S)]\ < |/l'V(jnfJ^2,iV+A3,„,w + A4,„,w + A5,„,Ar} + ^mf JA6,„,J+ (4.11) 



(ii) For any Lipschitz function h, and for every n > K , 

2 



E[/l(5„)] -E[/l(5)]| < l/l'loo ^^nf ^{^6.n,J+A7.,/} 



N>,n\\2as ((2X + i.™)E[/2(X)])^/'^ (Ef=™^l,„)^ ^' 



\nSn<z)]-¥{S <z)\ < —(— inf |A6,„,j + A7,j| 



( Hi ) For any z g M, and for every n > K , 

2/2 

N>,n\\2as ((2X + ^™)E[/2(X)])'/'^ iHt^^ln)' 

Proof of Proposition^ Let us introduce a similar notation to NPP. Consider the Hilbert space Jo — M"'^ with 
elements u = (ut,;,l < t < n,l < I < v) G and the scalar product {utj,ut'j')f, 

Tn''' (f, t'). The ^— fold tensor product and the symmetrized tensor product of arc denoted by S)®^ and Jo®^, 
respectively. Let L^(X„) denote the space of r.v.'s subordinated to the Gaussian vector X„ := (X„(t))i<t<„. 
Any element ^ G (X„) admits a chaotic expansion ^ — X)fco ^(^) ) ' '^hsre g(i) S Jo*^^ and the linear map- 
ping : ^ L2(X„) satisfies /(f) (g) = /(f)(syni(g)), Ej2)(g) = €!||sym(5)||^»f , and E[/(f)(5)/(f,)(5')] = 
0, £ 7^ 5(£) € i^*^^, (7(f/) S Jo**^ , where sym denotes the symmetrization operator. In particular, for any 
t = 1, . . . , n, fee Z"^, |fc| =: I we have Hj^{Xn{t)) — L(e^ (^ge{k)), where 

giik) sym(M^f" (gj-.-iguf^^*"') = ^ 6(v; fe) Mt,,;i <» • • • «> 
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and where b{v;k) = sym[b{v; k)] is the symmetrization of the function {l,...,^}^ 3 v = {vi, . . . ,V£) i— > 
b{v; k) := nj^=i l{v^ ^ r,ki + . . . + K-i < i < h + . . . + K). Thus, 5„ = n^i/^ ^"^^ j (X„(i)) admits the 
chaotic expansion 



oo _j n 

Sn^Y^hiM) With g'l:=—Y^ ^^(^) 
l=m t=l V£{l,...,v}t 



Ut,ve, 



where be{v) := E|A;|=£('^/(^)/^0&(«; fc) depend only on / e L2(X„(t)) = h^{X) and satisfy E/(2)(X) = 
^■X]t;e{i ly}*^ ^K''^)' i'^ NPP. It is important that here the g^'s are symmetric since the 6f(t;)'s are 
symmetric. Therefore E/^^^^(5") = Next, for > to consider the truncated expansion 



N 



Sn,N ■= Y, hf)(9£ 



Note that 



N 



N 



t=m t,t' = lt),Vg{l,...,,y}f 

^ TV n £ 

^E^! E E h{v)b,{v')J\r':-<\t^t'). 



e=m t,t' = l i=] 

Using |ri^'-''^(t,t')| < 6l(t - t') similarly as in NPP we obtain 

3 



E[/i(5„)] -E[/i(^".^)] < 2(2^ + ^"^)I'^"U(e[/2(X)] ^ E[/f,)(X)]j < -|/i"|oo^2.Ar. (4.14) 



For N > m, let Z„ jv be a centered Gaussian random variable with variance E5'^ ^ = X^fcm '''^n' with ct^ ,^ 
defined in (|4.3p . (Note that the last variance is slightly different from the variance of Zjv in (NPP, sec. 4.2).) 
Let D denote the Malliavin derivative in L2(X„), see NPP. Using ^-iE|li:)/(f)(.gP)|||, = ^!|1.9"|||»« 
dmH), as in (NPP, (4.46)) we obtain 



""^ -^l.n. see 



E[/i(Z„,jv)] -E[/i(5„,jv)] 



Next, using (NPP, (3.39)) 



N 



< 



,n.N + A^n.N ,n,N)- 



(4.15) 



< - \h" 
- 2 ' 



N 

E 



1 ^ 



To estimate the difference af^ — cr|, we use an interpolation identity from Houdre et al. (1998). Let 
(Xi, X2), {Wi,W2) be two (2z^)— dimensional Gaussian vectors with zero means and respective covari- 
ance matrices ^X,XJ] = E[W.,WJ] = /, i = 1,2, E[XiXl] = Ei, E[PFiVF|] = Eo. For a e [0,1] 
let {Xia,X2a) denote the "interpolated" Gaussian vector with zero mean and E[XiaXj^] — I,i = 
1, 2, E[XiaXl^] — (1— a)So+aSi. Let / £ L^(X) be a real function satisfying the conditions of Proposition 
H Then from (1.1), (1.3)) we obtain 



|Cov(/ (Xi) , / {X2) ) - Cov{f{W,),f{W2)) I 



J E[a/(Xi„)T(I]i-I]o)9/(J£:2a)] da 
< |9/|Ll|Si-So||, 



(4.16) 
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where df = {df /dx(^\ . . . ,df /Ox^'^^y G W. Let F„(r) := ^J',^, Cov(/ (X„(M)) , / (^n(i')) ) , r G [0,1] 
so that cr^ = F„(r)dr. Using (14.16^ . for 1 < J < n we can write jo"^ — o-g| < Ri{n, J) + R2{n, J), where 



Ri{n,J) 



I E |Cov(/(X„(M)),/(X„(M +.?•)) -Cov(/(W^.(0)),/(W-,(j)))|dr < 2A,^n,j, 



i?2(n,J) < 2E[/2(X)]j.'" E ^'"(^) = 2A7,j. 

|fc|>,/ 

We also have \al - EfcmCrLI = Efcjv+i c^l.n < 5^2,jv, as in gH. Therefore, Efcm cr£,„ - 4 
,n,J + 2^7, J + i A2^N, implying 



< 



E[/^(Z„,w)] -E[/i(5)]| < \h"\oo[Ae,n,j + A^j + ^A2,N) for 1<J< 



(4.17) 



Finally combining (|4.14p . ( (|4.15p . and (|4.17p results in 



[/l(^„)] -E[/l(5)]| < |/l"|oo(v42,iV+A3 



El . 

< |^"|oo( inf {A2N+A3n.N+A4n.N+A5n.N }+ inf {Ae,,,j + AjA 
\N>m - J l<J<n i- ' ■ ' 

proving the bound in (14. lip . 

(ii) Following (NPP, proof of Theorem 2.2-(2)) and the previous results, for a Lipschitz function h we obtain: 
\E[h{S,,)]~E[h{Sn,N)]\ < \h'\^{{2K + ,y'^)E[f{X)]y'^^A2,N, 

N 

,n,N + Ai,n,N + ^5,n,Ar j 

l=m 

and \E[h{ZnM)]-E[h{S)]\ < 1^ (1 A2,Ar + ^ mf JA6,„,j + Ay,,/)) 
and therefore (|4.12p is established. 



(iii) Bound (|4.13p is obtained exactly as in (NPP, proof of Theorem 2.2-(3)). 



□ 



5 Applications of Lemma [T] and Theorem [T] 
5.1 Application to the IR statistic 

This application was developed in Bardet and Surgailis (2011, 2012). Let (Xf)jg[Q,2] be a continuous time 
Gaussian process with zero mean and generally nonstationary increments locally resembling a fractional 
Brownian motion with Hurst parameter H{t) e (0, 1). Consider the Increment Ratio (IR) statistic 

n— 3 I A 2,n -r^ . \2,n 



and 



with Af,'"X = ^(fc+2)/n ~ 2X(fc_|_i)/„ + Xf^/n and the convention ^ :— 1. Let ct| „(fc) := E 

I72,«(fc) Cr2,„(fc) 

Then i?2^"(X) = ;r^Er=o /(^«(^)) ' /(2;('^a;('^) I-t^^) + xf^) |/(|a;(i) | + jx^^)]) can be written as the 
sum of nonlinear function / of Gaussian vectors Yn{k) — {Yy\k),Yn'^\k)) e M^, < fc < n — 3. These 
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Gaussian vectors can be standardized, leading to the expression — -;^z^'Y^X=X fn.k i-^nik)) of the 

IR statistics as the sum of some functions fn,k of standardized Gaussian vectors JC„(fc), < fc < n — 3. (It is 
easy to check that the centered functions fn.k — E[/n,fc(^)] have the Hermite rank 2.) If (Xt) satisfies some 
additional conditions (specifying the decay rate of correlations of increments and the convergence rate to the 
tangent process), Theorem [1] can be applied to establish that y^{R^^"{X) - A{H{t)) dt) Af{0, a^) 

n— >-oo 

with an explicit function A and a variance cr^. An application of Lemma [T] to bound the 4th moment 
E(i?^'"(X) — Ei?^'"'(X))^ provides a crucial step in the proof of the almost sure consistency of the IR 
statistic, i.e. i?2."(X) ^ J^A{H{t))dt. See Bardet and Surgailis (2011) for details. Local versions of 

n— >-oo 

the IR statistic for point- wise estimation of H{t) are developed in Bardet and Surgailis (2012). The study 
of the asymptotic properties of these estimators in the last paper is also based on Theorem [1] and Lemma [1] 

5.2 A central limit theorem for functions of locally stationary Gaussian pro- 
cesses 

Using an adaptation of Dahlhaus and Polonik (2006, 2009), we will say that {Xt^n)i<t<n,n&N* is a locally 
stationary Gaussian process if 

Xt,n ■■= at,n(j) £t-j, for aU 1 < t < n, n e N*, (5.1) 

where {ek)k£Z is a sequence of independent standardized Gaussian variables and for 1 < t < n ^ N* the 
sequences {at^n{j))jez are such that there exist K > and a < 1/2 satisfying for all n e N* and j E Z, 

K 

max \at^nij)\ < — , with Uj := max(l, |j j"^-"^) for j G Z (5.2) 

l<t<n ' Uj 

and such that there exist functions r € (0, 1] ^{Tjj) E R satisfying the following conditions: 

sup \a{T,j)\ < — , VjeZ, (5.3) 

tG(0,1] 

and sup max |(afe,„(j) - a(r, j)| -)■ 0, W j E Z, VL > 0. (5.4) 

re(0,l] \[-n-T]-k\<L 

For T E (0,1] introduce a stationary Gaussian process 

Writ):=J2a{r,j)et-j, t E Z. 

with spectral density griv) — |a(r, u)p, a(r, w) (27r)~^/^ ^^.^^ e~'''''a(T, j), v E [— tt, tt]. Let 
Y„(fc) {Xk+i,n,---,Xk+.,„y, WrU) := {WAj + 1), ■ ■ ■ ,Wrij + ly))' . 
Note (Wr(j))jez is a M''— valued stationary Gaussian process. Let 

:= E[Yn{k)Yn{ky], J:^ ■■^E[Wr{0)Wr{Oy]. 
Proposition 3 In addition to i5.1\) - ^5.4^ , assume that 

sup IIE^^II < oo. (5.5) 

1-5(0,1] 

Let /fc „ E LQ(y„(fc)), 1 < k < n,n > 1 be a triangular array of functions all having a generalized Hermite 
rank at least m > 1/(1 — 2a). Let there exists a lJ^{X) — valued continuous function (pT, t G (0, 1] such that 
relation HIM holds, with fuA^) fkA^]!^^)- Then the CLT of (E23j holds, with 

a^:=f dTY,H^r{Wr{0))(Pr{Wr{j))] (5.6) 

and 4>t{x) :— (firij^r x) defined as in Corollary\^ 
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Proof. We apply Corollary [21 Let us first check 



sup ||S[„^]_„ - E^ll — > 0. (5.7) 

re(Oa] n-s-oo 



We have 



|o-[nr],«(P,g) - Cr^(P,'7)| = \^{a[nT]+pAP + j)'^[nT]+qA<} + j)-"'i'^:P + j)"-i'^:1+j)) <Tn,J + T^^j, 

where 

TL,j — 2K'^ ^ Up+jUq+j, T/Zj := ^ |a[„^]+p_„(p + j)a[„^]+q^„((7 + j)-a(r,p + j)a(T,q+j))| 

m>j \j\<'j 

according to (15.21) and (15.31) . Clearly, T,' j can be made arbitrarily small by choosing J large enough. Then for 
any J < oo fixed, we have that sup^^j-Q ^ T"" j according to assumption (j5.4p . This proves (|5.7p . In a simi- 
lar way, one verify that for any r e (0,1], j,/ G Z, ||E[Y-„([nT]+j)r„([7iT]+j')^]-E[W^,(j)T^,(i')^] II ^ 
implying condition (13. 8p . The dominating condition (|3.3I) on cross-covariances is ensured by (|5.2p and the 
fact that (1 — 2a)r7i > 1. The remaining conditions of Corollary [2] are trivially satisfied. □ 



Remark 2 Dahlhaus and Polonik (2006, 2009) discussed the short-memory case {at,nU))j&^ G 1 < t < n 
only. On the other hand, condition (|5.2p allows for the long-memory case {at^n{j))j& G ht,«(j)l = 

oo. The last case is also discussed in Roueff and von Sachs (2010), where similar conditions as (|5.2p and 
(|5.3p are provided in spectral terms. It is not clear whether condition (j5.4p allows for jumps of the parameter 
curves r i— > a(T, •) as in Dahlhaus and Polonik (2006, 2009), in particular, for abrupt changes of the memory 
intensity of Gaussian process (|5.ip . See also Lavancier et al. (2011) for a related class of nonstationary 
moving average processes with long memory. 

Remarks Note that ccTE^a; = /_^^ .g^(w)| ^J^^ e'^'^x^^^ |^dw for any cc = (.t^^^ . . . , 2:(''))t e R''. Therefore 
condition 'vLdve[--K,iT],Te(o,i]9r{v) > 7 > on the spectral density of {Wr{t)) implies condition (j5.5p . since 
x'^Y^rX > c|a;p, c := 27rz^7 > 0. 

Remark 4 For stationary Gaussian long memory process, condition mil — 2a) > 1 was first obtained in 
Taqqu (1975). Proposition [3] can be applied to prove the asymptotic normality of various statistics of locally 
stationary processes, see, e.g., Roueff and von Sachs (2010). 

Acknowledgment. The authors are grateful to two anonymous referees for valuable suggestions and com- 
ments that helped to improve the original version of the paper. 
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